Irreversible Magnetization of Pin-Free Type II Superconductors 
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The magnetization curve of a type II superconductor in general is hysteretic even when the 
vortices exhibit no volume or surface pinning. This geometric irreversibility, caused by an edge 
barrier for flux penetration, is absent only when the superconductor has precisely ellipsoidal shape 
or is a wedge with a sharp edge where the flux lines can penetrate. A quantitative theory of 
this irreversibility is presented for pin-free disks and strips with constant thickness. The resulting 
magnetization loops are compared with the reversible magnetization curves of ideal ellipsoids. 



PACS numbers; 74.60. Ec, 74.60.Ge, 74.55.+h 

The magnetic moment of most superconductors is well 
known to be irreversible. After Abrikosov's ||l| prediction 
of quantized flux lines it became clear |^ that the mag- 
netic hysteresis is caused by pinning of these vortex lines 
at inhomogeneities in the material. Flux-line pinning and 
the related critical state |^ were subsequently confirmed 
quantitatively in numerous papers [Q. However, similar 
hysteresis effects were also observed |^ in type I super- 
conductors, which do not contain flux lines but normal 
conducting domains, and in type 11 superconductors with 
negligible pinning. In these two cases the magnetic ir- 
reversibility is caused by a geometric (specimen-shape 
dependent) barrier which delays the penetration of mag- 
netic flux but not its exit. In this respect the geometric 
barrier behaves similar to the Bean-Livingston barrier 
[^,0 for vortices penetrating a parallel surface. 

The geometric irreversibility is most pronounced for 
thin films of constant thickness in a perpendicular field. 
It is absent only when the superconductor is of exactly 
ellipsoidal shape or is tapered like a wedge with a sharp 
edge where flux penetration is facilitated. In ellipsoids 
the inward directed driving force exerted on the vortex 
ends by the surface screening currents is exactly compen- 
sated by the vortex line tension |^ , and thus the magneti- 
zation is reversible. In specimens with constant thickness 
(i.e. rectangular cross-section) this line tension opposes 
the penetration of flux lines at the four corner lines, thus 
causing an edge barrier; but as soon as two penetrating 
vortex segments join at the equator they contract and are 
driven to the specimen center by the surface currents, see 
Fig. 1 below. As opposed to this, when the specimen pro- 
file is tapered and has a sharp edge, the driving force even 
in very weak applied field exceeds the restoring force of 
the line tension such that there is no edge barrier. The 
resulting absence of hysteresis in wedge-shaped samples 
was nicely shown by Morozov et al. . 

An elegant analytical theory of the field and current 
profiles in thin superconductor strips with an edge bar- 
rier has been presented by Zeldov et al. see also the 
extensions ||ll|]. With increasing applied field Ha, the 
magnetic flux does not penetrate until an entry field H^n 



is reached; at Ha = H^n the flux immediately jumps to 
the center, from where it gradually fills the entire strip 
or disk. This behavior in increasing Ha is similar to that 
of thin films with artificially enhanced pinning near the 
edge 1 1 1 1^ , but in decreasing Ha the behavior is differ- 



ent: In films with enhanced edge pinning (critical current 
density Jc,cdgo) the current density J at the edge imme- 
diately jumps from + Jc,cdgc to —Jc, edge when the ramp 
rate inverses sign, while in pin-free films with geometric 
barrier the current density at the edge first stays con- 
stant or even increases and then gradually decreases and 
reaches zero at Ha = 0. The entry field Hen was esti- 
mated for pin-free thin strips in Refs. pO|,p^, see also 
Refs. 
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In this letter the geometry-caused magnetic irre- 
versibility of ideal pin-free type II superconductors is cal- 
culated and discussed for the two most important exam- 
ples of circular disks (or cylinders) and long strips (or 
slabs) with rectangular profile of arbitrary aspect ratio 
b/a. 1 present flux-density proflles and magnetization 
loops and give explicit expressions for the entry held H^n 
and for the reversibility field -ffjcv above which the mag- 
netization curve is reversible. Finally, the modification 
of these results by volume pinning is briefiy mentioned. 

Let us first consider the magnetization of ideal el- 
lipsoids. If the superconductor is homogeneous and 
isotropic, the magnetization curves M{Ha', N) are re- 
versible and may be characterized by a demagnetizing 
factor N with < < 1. If Ha is along one of the three 
principal axes of the ellipsoid then is a scalar. One 
has A^ = for long specimens in parallel held, A^ = 1 
for thin films in perpendicular field, and A^ = 1/3 for 
spheres. If the magnetization curve in parallel field is 
known, M{Ha; 0) = B/ fio — Ha where B is the flux den- 
sity or induction inside the ellipsoid, then the homoge- 
neous magnetization of the general ellipsoid, M{Ha;N), 
follows from the implicit equation 



H,^Ha~NM{Hf,0). 



(1) 



Solving Eq. (1) for the effective internal held Hi, one ob- 
tains M = M{Ha;N) = M{Hi;0). In particular, for the 
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Meissner state {B = 0) one finds M{Ha; 0) = —Ha and 

Ha 



M(Ha;N) = - 



1 - N 



for \Ha\ < {l-N)H,i. (2) 



At the lower critical field Hd one has Hi = Hd, Ha = 
H'd = (1 - N)Hd, B ^ 0, and M = -Hd- Near the 
upper critical field Hc2 one has an approximately linear 
M{Ha] 0) = 7(7?a - i?c2) < with 7 > 0, yielding 



M{Ha;N) 



7 



l + 'jN 



{Ha-H,2) for Ha^H,2. (3) 



Thus, if the slope 7 ^ 1 is small (and in general, if 
\M/Ha\ ^ 1 is small), demagnetization effects may be 
disregarded and one has M{Ha \ N) k, M{Ha \ 0). 

The ideal magnetization curve of type II superconduc- 
tors with iV = 0, M{Ha] 0) or B{Ha; 0) = Ha+M{Ha] 0), 
may be calculated from Ginzburg-Landau (GL) theory 
[^6|, but any other model curve may be used provided 
M{Ha;0) — —M{—Ha\0) has a vertical slope at Ha — 
Hci and decreases monotonically in size for Ha > Hd- 
For simplicity in this letter I shall assume Hd «C Hc2 
(i.e. large GL parameter k ^ 1) and Ha <C Hc2- To il- 
lustrate the essential features I may thus use the realistic 
model M{Ha;0) = -Ha for \Ha\ < Hd and 
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Pin-free superconducting slab and strip 



FIG. 1. The magnetic field lines of B{x,y) in slabs or 
strips with aspect ratio b/a = 2 (top) and b/a — 0.3 
(bottom) in perpendicular magnetic field Ha- Top left: 
Ha/Hci — 0.66, in increasing field shortly below the entry 
field Hcn/Hci = 0.665. Top right: Ha/Hci = 0.5, decreasing 
field. Bottom: Ha/ Hd — 0.34 in increasing field just above 
Hen/Hci = 0.32. The field lines of cylinders look very similar. 
Note the straight field lines in the corners, corresponding to 
flux lines under tension. 



M{Ha;0) = {Ha/\Ha\){\Ha\'' - Hl^f'' - Ha (4) 

for \Ha\ > Hd, see the curve labeled oo in Fig. 3 below. 

In nonellipsoidal superconductors the induction B(r) 
in general is not homogeneous, and so the concept of 
a demagnetizing factor does not work. However, when 
the magnetic moment m = ^ J r x J(r)d"^r is directed 
along Ha, one may define an effective demagnetizing fac- 
tor N which in the Meissner state {B = 0) yields the 
same slope M/Ha = -1/(1 - N), Eq. (2), as an el- 
lipsoid with the same volume V. Here the definition 
M — m/V with m = mHa/Ha is used. For long 
strips and circular disks or cylinders with cross-section 
2a X 26 in a perpendicular or axial magnetic field along 
the thickness 2b, approximate expressions for the slopes 
M/Ha = m/{VHa) are given in Refs. [|7|jl|]. Using this 
and defining q = {\AI/Ha\ — l){b/a), one obtains the ef- 
fective iV for any aspect ratio b/a in the form 
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In the limits b <^ a and b ^ a, formulae (5) are exact, 
and for general b/a the relative error is < 1%. For a — b 
(square cross-section) they yield for the strip N — 0.538 
(while N = 1/2 for a circular cylinder in perpendicu- 
lar field) and for the short cylinder N — 0.365 (while 
iV 1/3 for the sphere). 

Next we consider the full, irreversible magnetization 
curves M{Ha) of pin-free strips and cylinders with cross 
section 2a x 2b. Appropriate continuum equations and 
algorithms (which apply also to pinning) have been pro- 
posed recently by Labusch and Doyle |^9| and by the 
author , based on the Maxwell equations and on con- 
stitutive laws which describe fiux fiow and pinning [or 
thermal depinning expressed, e.g., by an electric field 
E(J, B)] and the reversible magnetization in absence of 
pinning, AI{Ha\ 0). Here I shall use the method [20] and 
the model M{Ha;0), Eq. (4). The pin-free flux dynam- 
ics will be described as viscous motion by E = ppf{B)J 
with flux-flow resistivity ppF oc _B. In both methods the 
M {Ha; 0) law enters the driving force density on the vor- 
tices, Jh X B with definition Jh = V x H, where H(B) is 
obtained by inverting the relation B(H) = H-|-M(H; 0). 

While method [19] considers a magnetic charge den- 
sity on the specimen surface which causes an effective 
field Hi(r) inside the superconductor, our method [20] 
couples the arbitrarily shaped superconductor to the ex- 
ternal field B(r,t) via surface screening currents: In a 
first step the vector potential A(r, t) is calculated for 
given current density J; then this relation (a matrix) is 
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FIG. 2. Top: Profiles of tfie axial magnetic induction 
By{r,y) in the midplane y = of a pin-free superconductor 
disk with aspect ratio b/a — 0.3 in increasing field (solid lines) 
and then decreasing field (dashed lines), plotted at Ha/Hci 
= 0.4, 0.42, . . ., 0.5, 0.52, 0.6, 0.5, 0.7, 0.8, 0.7, 0.6, . . ., 0.1, 0. 
Bel = i-ioHci- Bottom: The induction By (0,0) in the center 
of the same disk (solid line) and of a strip (dashed line), both 
with b/a = 0.3. The symbols mark the field values at which 
the profiles are taken. Also shown is the magnetization loop 
for the same disk and strip and the corresponding reversible 
magnetization (dotted lines), see also Fig. 3. 

inverted to obtain J for given A and given Hq; next the 
induction law is used to obtain the electric field [in our 
symmetric geometry one has E(J,B) = —dA/dt], and 
finally the constitutive law E = E(J, B) is used to elimi- 
nate A and E and obtain one single integral equation for 
J(r, t) as a function of Ha(i), without having to compute 
B(r,i) outside the specimen. This method in general is 
fast and elegant; but so far the algorithm is restricted 
to moderate aspect ratios, 0.03 < b/a < 30, and to a 
number of grid points not exceeding 1000 (on a Personal 
Computer) . Improved accuracy is expected by combining 
methods (19) (working best for small b/a) and (20). 

The penetration and exit of flux computed by method 
[20] is illustrated in Figs. 1 and 2 for isotropic strips and 
disks without volume pinning, using a flux-flow resistiv- 
ity pff = pB{r) with p ~ 140 (strip) or p = 70 (disk) in 
units where Hd = a = p,Q = \dHa/dt\ ~ 1. The profiles 
of the induction By (r, y) taken along the midplane y — Q 
of the thick disk in Fig. 2 have a pronounced minimum 
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FIG. 3. Irreversible magnetization curves —M{Ha) of 
pin-free circular disks or cylinders with aspect ratios b/a = 
0.08, 0.15, 0.25, 0.5, 1, 2, 5, and oo in axial field (solid lines). 
In these type II superconductors the irreversibility is due to 
a purely geometric edge barrier for flux penetration. The 
dashed curves are the reversible magnetization curves of the 
corresponding ellipsoid defined by Eqs. (1,4,5). 

near the edge r — a, precisely in the region where strong 
screening currents flow. Away from the edges, the current 
density J = V x B//io is nearly zero; note the parallel 
field fines in Fig. 1. The quantity Jh = V x H(B) which 
enters the Lorentz force density Jh x B, is even exactly 
zero since we assume absence of pinning. Our finite flux- 
flow parameter p and finite ramp rate dHa/dt = ±1 mean 
a dragging force which, similar to pinning, causes a weak 
hysteresis and a small remanent flux at Ha = 0; this ef- 
fect may be reduced by choosing larger resistivity and 
slower ramping. 

The induction _By(0, 0) in the specimen center in Fig. 
2 performs a hysteresis loop very similar to the magneti- 
zation loops M{Ha) shown in Figs. 2, 3. Both loops are 
symmetric, e.g., M{—Ha) = —M{Ha)- The maximum of 
M{Ha) defines a field of first flux entry i/cn, which closely 
coincides with the field at which i3j,(0, 0) starts to 
appear. The computed entry fields are well fitted by 



Hf/'P/H,i = tanh ^0.366/0, 

H^^^/Hci = tanh ^/OMbJ^ . (6) 

These formulae are good approximations for all aspect ra- 
tios < 6/a < cx), see also the estimates of i/on ~ \Jb/a 
for thin strips in Rcfs. [|lO|,0. 

The virgin curve of the irreversible M{Ha) of strips 
and disks at small Ha coincides with the ideal Meissner 
straight line M = —Ha/ (1 — N) of the corresponding el- 
lipsoid, Eqs. (2,5). When the increasing Ha approaches 
Hem flux starts to penetrate into the corners in form of 
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FIG. 4. Magnetization curves of a thicli disk with aspect 
ratio hja = 0.25 for various degrees of volume pinning, = 
0, 0.25, 0.5, 1, 1.5, 2, 3, 4 in units Hci/a, and for various 
sweep amplitudes. The inner loop belongs to the pin-free 
disk (Jc = 0), the outer loop to strongest pinning. Also 
shown is the reversible magnetization curve of the corre- 
sponding ellipsoid (dashed curve). All loops are symmetric, 
M(-H„) = ~M{Ha). 

stretched flux lines (Fig. 1) and thus \M{Ha)\ faUs be- 
low the Meissner line. At Ha = Hen flux penetrates and 
jumps to the center, and \M{Ha)\ starts to decrease. In 
decreasing Ha, this barrier is absent. As can be seen 
in Fig. 3, above some field i?rcv, the magnetization curve 
M{Ha) becomes reversible and exactly coincides with the 
curve of the ellipsoid defined by Eqs. (1, 4, 5) (in the qua- 
sistatic limit with p~^dHa/dt 0). The irreversibility 
field H^cv is difficult to compute since, in our present 
algorithm, it slightly depends on the choices of the flux- 
fiow parameter p (or ramp rate) and of the numerical 
grid, and also on the model for M{Ha\ 0). In the interval 
0.08 < h/a < 5 wc find with relative error of 3%, 

H^^^/Hci = 0.65 + 0.12 In (b/a) , 

H^l^/H,i ^ 0.75 + 0.151n(6/a) ■ (7) 

This fit obviously does not apply to b/a ^ 1 (since -ffrov 
should exceed H^n > 0) nor to 5/a ^ 1 (where i/rov 
should be close to Hd). The limiting value of TJrov for 
thin films with 6 <C a is thus not known at present. 

Remarkably, the irreversible magnetization curves 
M{Ha) of pin- free strips and disks fall on top of each 
other if the strip is chosen twice as thick as the disk, 
(&/a)strip ~ 2(6/a)disk- This striking coincidence holds 
for all aspect ratios < 5/a < oo and can be seen 
from each of Eqs. (5-7): The effective N [or virgin slope 
1/(1 — N)], the entry field iJcn, and the reversibility field 
i/i-cv are nearly equal for strips and disks with half thick- 
ness, or for slabs and cylinders with half length. 

Another interesting feature of the pin-free magnetiza- 
tion loops is that the maximum of \M{Ha)\ exceeds the 



maximum of the reversible curve (equal to Hd) when 
b/a < 0.8 for strips and b/a < 0.4 for disks, but at larger 
b/a it falls below Hd- The maximum magnetization may 
be estimated from the slope of the virgin curve 1/(1 — A^), 
Eq. (5), and from the field of first flux entry, Eq. (6). 

Finally, Fig. 4 shows how the irreversible magnetiza- 
tion loop is modifled when volume pinning of the flux 
lines is switched on. Increasing critical current density 
Jc (in natural units Hd / a) inflates the loops nearly sym- 
metrically about the pin-free loop or (above i?rev) about 
the reversible curve, and the maximum of \M{Ha) \ shifts 
to higher fields. Above Hrov the width of the loop is 
nearly proportional to Jc, as expected from previous the- 
ories QJl^,|l^ which assumed Hd = 0, but at small fields 
the influence of finite Hd is clearly seen up to rather 
strong pinning. 
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